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SYNOPSIS 


The  earth's  upper  atmosphere  continuously  emits  infrared  radiation. 
Some  of  this  radiation  is  powered  by  the  energy  stored  in  vibrationally 
excited  nitrogen  that  can  be  produced  by  a  number  of  exothermic  pro¬ 
cesses.  One  of  these  processes  is  the  reaction  NO  +  N  -*  N2  +  0,  in  which 
3.27  eV  is  available  for  vibrational  and  rotational  energy  of  the  N2 
product  as  well  as  for  kinetic  energy. 

The  first  task  of  the  research  program  was  to  determine  the  fraction 
of  the  energy  appearing  in  Na  vibrations  by  the  use  of  Raman  spectroscopy 
to  determine  the  population  of  Na(v=l)  that  is  produced.  These  measure¬ 
ments  established  that  2S  ±  3$  of  the  exothermicity  of  the  reaction  ap¬ 
pears  as  vibrational  energy  in  the  nitrogen  and  that  the  NO  +  N  -*  Na  +  0 
reaction  is  an  important  source  of  vibrationally  excited  Na  in  the  E 
region.  An  article  describing  the  measurements  has  been  accepted  by  the 
Journal  of  Chemical  Physics  and  is  attached  as  Appendix  A  to  this  report. 

Once  produced,  the  vibrationally  excited  nitrogen  is  deexcited  only 
by  collisions.  In  the  E  region,  one  important  quenching  mechanism  is 
energy  transfer  to  C0a,  the  rate  for  which  is  moderately  well  known.  A 
second  important  deexcitation  mechanism  was  considered  to  be  V-T  energy 
exchange  with  0(*P).  At  the  start  of  this  research  program,  the  rate 
coefficients  for  this  latter  process  were  only  known  for  translational 
temperatures  far  higher  than  are  found  in  the  earth's  atmosphere. 

Therefore,  Task  II  of  our  work  was  to  measure  this  rate  down  to  as 
low  a  temperature  as  possible  by  the  use  of  a  shock  tube  with  ozone  as 
the  source  of  atomic  oxygen  and  the  infrared-tracer  method  to  monitor 
the  vibrational  energy  in  the  nitrogen.'  These  measurements  have  estab¬ 
lished  that  V-T  relaxation  of  Na{v=l)  by  0(3P)  is  an  efficient  process 
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with  a  rate  varying  from  1  atm-y,sec  at  3000  K  to  2  atm-^sec  at  1200  K. 
Our  results  are  in  agreement  with  recent  results  at  lower  temperature 
in  an  afterglow  experiment  carried  out  by  McNeal,  Whitson,  and  Cook.1 
In  addition,  the  results  show  that  0(3P)  V-T  relaxation  plays  an  im¬ 
portant  role  in  determining  the  vibrational  temperature  of  Na  in  the  E 
region.  The  results  are  fully  described  in  an  article  that  has  been 
submitted  for  presentation  and  publication  in  the  proceedings  of  the 
Ninth  International  Shock  Tuba  Symposium  (Stanford  University,  July 
1973).  This  article  constitutes  Appendix  B  to  this  report. 

In  summary,  quenching  of  vibrationally  excited  nitrogen  by  atomic 
oxygen  appears  to  be  an  important  loss  process  for  Na(v=l)  in  the  lower 
thermosphere.  The  use  of  correct  rates  for  this  process  results  in  a 
large  reduction  in  the  expected  Na  vibrational  temperature  in  this 
region. 
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I  INTRODUCTION 


In  the  E  region  of  the  earth's  atmosphere,  a  number  of  exothermic 
chemical  reactions  occur  that  can  produce  Na,  vibrationally  excited 
nitrogen.  In  the  lower  E  region  (below  125  km),  a  major  loss  of  nitro¬ 
gen  vibrational  energy  is  by  energy  transfer  to  00a, 

Na ( v=l )  +  C0a  -  COa(v3=l)  +  Na  -  18cm"1  ,  (l) 

with  subsequent  radiation  izt  the  infrared  (principally  at  4.3  p.) •  Pro¬ 
cesses  that  excite  or  deexcite  Na  are  thus  of  particular  interest  in 
connection  with  some  defense  problems  concerned  with  atmospheric  infra¬ 
red  measurements. 

* 

The  main  processes  that  have  been  considered  for  producing  Na  are: 


N(4S)  +  NO  -  Na  +  0(3P)  +  3,27  eV  (2) 

OOD)  +  Na  -  Na  +  0(3P)  +  1.96  eV  (3) 

Na(A3Z*)  +  0(3P)  -  Na  +  0(3P,*D,‘S),  AE  =  6.2,  4.2,  2.0  eV  (4) 

Na(A3l£)  +  0a  -  Na  +  0a(9Z,1 AE  ■  6.2,  5.2,  4.6  eV  (5) 

e"  +  Na  -  (Na)  -  Na  +  e"  .  (6) 


Also,  emission  of  Vegard-Kaplan  band  radiation  by  Na(A3E*)  in  most 
instances  will  leave  Na  vibrationally  excited.  Whereas  process  (6)  is 
cortainly  important  in  auroral  conditions,  it  is  not  important  in  the 
undisturbed  E  region. 

If  all  these  processes  wore  equally  effective  (as  a  fraction  of 
the  energy  available)  in  producing  Na,  then,  in  the  E  region,  process 


1 


(3)  would  produce  more  N2  than  process  (2),  and  the  latter  would  produce 
more  than  the  sum  of  processes  (4)  and  (5).  The  first  goal  of  the  studies 
at  SRI,  therefore,  has  been  to  investigate  processes  (2)  and  (3). 

A  second  aspect  of  the  atmospheric  processes  powering  emissions 
$ 

through  N2  concerns  the  reactions  that  can  quench  it  in  the  upper  atmo¬ 
sphere.  The  most  important  deexcitation  mechanisms  have  been  considered 
to  be  energy  transfer  to  C02  in  the  lower  E  region  [process  (l)] — the 
rate  for  which  is  moderately  well  known — and  diffusion,  electron  quench¬ 
ing,  and  V-T  energy  exchange  with  0(3P)  in  the  upper  E  region  and  higher. 
At  the  time  this  program  was  undertaken,  rate  coefficients  for  the  latter 
process  had  been  measured  only  in  the  temperature  range  from  3000  to 
4500  K.2  A  second  goal  of  our  work  (Task  II )  has,  therefore,  been  mea¬ 
surement  by  shock  tube  techniques  of  this  rate  for  temperatures  down  to 
1200  K.  Recently,  low  temperature  measurements  of  the  rate  of  vibrational 
relaxation  of  Na  by  0(3P)  have  been  reported  which  are  many  orders  of  mag¬ 
nitude  faster  than  was  anticipated  on  the  basis  of  theoretical  extrapola¬ 
tions  of  the  high  temperature  results.1  The  intermediate  temperature 
results  of  this  study  are  therefore  of  importance  to  verify  the  low  tem¬ 
perature  data  and  to  establish  the  correct  temperature  dependence  of  this 
rate. 
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II  TASK  I:  VIBRATIONAL  POPULATION  MEASUREMENTS 


A.  Background 

One  aim  of  the  Task  I  program  has  been  to  determine  the  efficiency 
with  which  processes  (2)  and  (3)  produce  N2  . 

Measurements  of  process  (2),  the  N  +  NO  reaction,  a  decade  ago  indi¬ 
cated  that  about  30$  of  the  exothermicity  of  the  reaction  appears  in  N*  ?“5 
If  this  is  true,  and  if  the  efficiency  of  (3)  is  as  low  as  was  indicated 
by  the  theoretical  study  of  E.  R.  Fisher  and  E.  Bauer6,  then  process  (2) 
is  much  more  important  than  (3)  as  a  source  of  N2  .  For  these  reasons, 
we  proposed  that  a  measurement  of  the  efficiency  of  production  of  N2  in 
process  (2)  should  be  a  primary  part  of  the  current  program. 

We  had  also  intended  to  make  measurements  of  process  (2)  by  the 
infrared-tracer  technique  with  CO.  Unknown  to  us  at  the  time  of  writing 
the  proposal,  such  a  program  was  already  under  way  at  Air  Force  Cambridge 
Research  Laboratories.  That  program  revealed  serious  problems  in  inter¬ 
pretation  when  the  infrared-tracer  technique  was  applied  to  process  (2). 
Early  in  March  1972,  the  Optical  Physics  Group  of  A.  T.  Stair  at  AFCRL 
discussed  their  findings  with  us.  The  problems  have  also  been  discussed 
with  Drs.  Gilmore  and  Bauer.  In  all  these  discussions,  general  agreement 
has  been  reached  that  process  (2)  should  be  given  high  priority  and  that 

the  preferable  technique  for  studying  it  would  be  the  Raman  scattering 

* 

technique  for  measuring  N2  .  The  Raman  technique  was  anticipated  to  be 
easier  to  apply  to  process  (2)  than  it  was  in  our  earlier  work  on  pro¬ 
cess  (3), 7  both  because  much  higher  initial  populations  of  N  atoms  could 
be  created  than  was  possible  with  0(*D)  and  because  of  the  expectation 
that  about  30$  of  the  reaction  energy  would  appear  as  vibrations  in 
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Although  the  main  goal  of  the  measurements  was  to  obtain  reliable 
information  on  the  N  +  NO  reaction,  it  was  hoped  that  the  use  of  the 
Raman  technique  for  reaction  (2)  would  reveal  the  causes  of  the  diffi¬ 
culties  encountered  by  others  in  using  the  infrared-tracer  technique  on 
,that  reaction.  In  this  way,  it  should  be  possible  to  improve  the  reli¬ 
ability  of  the  infrared  technique. 

Last  year  our  measurements  were  confined  to  reaction  (3) — the 
quenching  of  O^D)  by  nitrogen.  The  vibrationally  excited  nitrogen  was 
measured  by  Raman  spectroscopy.  During  the  period  of  that  experimental 
program,  a  theoretical  study  by  E.  R.  Fisher  and  E.  Bauer  concluded  that 
less  than  5$  of  the  0(*D)  electronic  energy  is  channeled  into  the  vibra¬ 
tional  levels  of  Na  . 6  Our  experimental  findings  were  inconclusive  (see 
reference  10  of  Appendix  A)  and  made  further  work  on  this  reaction  de¬ 
sirable.7  One  objective  of  the  current  program  was  to  perform  such  mea- 

* 

surements  by  using  CO  as  an  infrared-emitting  tracer  of  the  N2  .  However, 
because  of  the  difficulties  in  interpreting  the  results  when  the  infrared- 
tracer  technique  was  applied  to  process  (2),  its  application  to  process  (3) 
was  postponed  to  some  future  time. 

B.  Accomplishments 
The  reaction 

N(4S)  +  NO  -  N2( v=n)  +  0(3P)  (2) 

is  fast  («2xl0“11  cm3  molec”1  sec”1  at  room  temperature8*9)  and  is  suf¬ 
ficiently  exothermic  (3.27  eV)  to  produce  N2  up  to  n=12.  It  is  expected 
that  the  initial  vibrational  distribution  will  be  rapidly  degraded  by 
vibrational  exchange  collisions  of  the  type 

N2(v=n)  +  N2(v=0)  -«N2(v=n')  +  N2(v=n-n')  ,  (7) 
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resulting  in  an  increased  population  of  N2(v=l) . 5-7  > 10“15  Hence  a  mea¬ 
surement  of  N2(v=l)  should  be  sufficient  to  determine  the  fraction  of 
the  exothermicity  of  reaction  (2)  that  appears  as  vibrational  energy. 

O 

Such  measurements  have  been  made  by  observing  Raman  scattering  of  4880  A 

Ar+  laser  radiation.  The  intensity  of  the  Q-branch  of  the  anti-Stokes 

1 

O 

transition  (at  4382  A)  was  measured,  because  the  anti-Stokes  spectra  are 
free  from  Raman  transitions  from  the  N2(v=0)  level  [whereas  the  Stokes 
transitions  from  N2(v=l)  lie  in  the  rotational  structure  of  the  Stokes 
lines  from  N2(v=0)]. 

The  measurements  are  fully  described  in  an  article  that  is  attached 
as  Appendix  A  to  this  report.  The  conclusion  reached  was  that  25  ±  3$ 
of  the  exothermicity  of  reaction  (2)  appears  as  vibrational  energy.  Of 
necessity,  the  use  of  a  rotational  Raman  spectrum  to  determine  gas  tem¬ 
perature  is  treated  only  briefly  in  Appendix  A.  The  purpose  of  the  next 
section  is  to  describe  this  measurement  in  more  detail. 

C.  Temperature  Determination  from  Rotational  Raman  Spectra 

With  heating  tapes  wrapped  around  the  pyrex  cell  (the  first  figure 

.  O 

of  Appendix  A),  nitrogen  at  atmospheric  pressure  was  heated  to  ~  100  C. 

Under  these  conditions,  with  a  thermocouple  close  to  the  focus  of  the 

0 

laser  beam  indicating  a  gas  temperature  of  100  C,  the  rotational  Raman 
spectrum  shown  in  Fig.  1  was  recorded.  The  series  of  lines  shifted  to 

O 

longer  wavelengths  of  the  4880  A  laser  line  are  called  the  Stokes  lines 
(j  -*  J  +2)  and  the  series  of  lines  shifted  to  shorter  wavelengths  are 
called  the  anti-Stokes  lines  (J  +  2  -*  J).  They  are  also  referrad  to  by 
some  authors  as  the  S  branch  and  0  branch,  respectively  [consistency 
with  international  nomenclature  would  require  that  both  series  of  lines 
be  called  S  branches,  since  AJ  =  j'-J"  (<J#  refers  to  the  upper  state  and 
J"  refers  to  the  lower  state)  is  +  2  in  both  cases].16  In  Fig.  1,  the 
numbering  of  the  lines  is  the  J  value  of  the  lower  state.  The  intensity 
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In  a  recently  published  paper, 18  gas  temperatures  were  deduced  by 
this  same  technique  from  observations  of  only  the  Stokes  lines  without 
correcting  for  the  v3  dependence  of  the  intensity.  The  tenperatures 
deduced  were  always  lower  than  those  measured  by  a  thermocouple.  This 
agrees  with  our  observations  if  the  v3  term  is  not  included.  This 
shows  the  importance  of  including  the  v3  variation  if  the  correct  tem¬ 
perature  is  to  be  obtained.  It  should  also  be  stressed  that  it  is  im¬ 
perative  to  determine  the  gas  temperature  from  both  branches  of  the 
rotational  Raman  spectrum  to  verify  the  constancy  of  the  counting  ef¬ 
ficiency  over  the  wavelength  range  of  the  spectrum. 
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Ill  TASK  II:  VIBRATIONAL  RELAXATION  OF  Na  BY  0(3P) 


A.  Background 


The  work  on  Task  II,  which  involves  the  measurement  in  a  shock  tube 
of  the  rate  of  vibrational  relaxation  of  Na  by  0(3P)  atoms,  is  a  continu¬ 
ation  of  work  carried  out  in  1970-71  under  ARPA  Contract  DAHC04-70-C- 
0036.  The  objective  has  been  to  determine  the  relaxation  time  p*rN  Q 
over  the  temperature  range  from  approximately  1,000  to  greater  than 
3,000  K  for  the  purpose  of  comparing  the  data  of  Breshears  and  Bird2  in 
the  range  from  3,000  to  4,500  K  and  establishing  the  temperature- 
dependence  of  the  rate  for  extrapolation  to  the  lower  temperatures  found 
in  the  upper  atmosphere. 


The  high-temperature  shock  tube  data  of  Breshears  and  Bird  was  ob¬ 
tained  using  the  laser-schlieren  technique,  in  which  the  density  gradi¬ 
ent  in  the  relaxation  region  behind  the  shock  wave  is  measured  through 
its  effect  on  the  refractive  index  of  the  gas.  It  was  originally  in¬ 
tended  that  the  same  technique  be  used  at  the  lower  temperatures  of 
this  study;  however,  it  was  found  that  at  the  lower  temperatures  the 
density  gradients  are  simply  too  small  and  disturbances  from  boundary 
layer  turbulence  are  too  large  for  the  schlieren  or  other  methods  rely¬ 
ing  on  refractive  index  changes  to  be  successful.  Therefore,  efforts 
were  concentrated  on  the  infrared  tracer  measurement  technique,  in  which 
the  radiation  from  small  amounts  of  CO  added  to  the  gas  mixture  is  pro¬ 
portioned  to  the  vibrational  energy  of  the  Na  . 


B.  Accompl. ghments 


The  experimental  program  thus  carried  out,  its  interpretation,  and 
the  results  are  described  in  full  in  Appendix  B,  which  is  a  preprint  of 


an  article  submitted  for  presentation  and  publication  in  the  proceed¬ 
ings  of  the  Ninth  International  Shock  Tubs  Symposium  to  be  held  at 
Stanford  University  in  July  of  1973.  As  described  in  Appendix  B,  the 
objective  of  Task  II  has  been  met  in  full  in  that  the  rate  of  vibra¬ 
tional  relaxation  of  Na  by  0(3P)  has  been  determined  for  the  tempera¬ 
ture  range  from  1200  to  3000  K.  In  addition,  some  limited  results  for 
the  rate  of  vibrational  relaxation  of  CO  by  0(3P)  in  the  same  tempera¬ 
ture  range  have  been  obtained. 


The  measured  values  of  pT^  Q  range  from  approximately  1  atm-y,sec 
at  3000  K  to  2  atm-p,sec  at  1200  K,  with  an  uncertainty  of  ±  50$.  Both 
the  magnitude  and  the  trend. of  the  results  are  in  excellent  agreement 
with  the  results  of  Breshears  and  Bird,2  as  well  as  with  the  recent  low 
temperature  measurements  of  McNeal,  et  al.1  In  marked  contrast  to 
theoretical  predictions,3  these  combined  results  establish  a  slowly 
varying  temperature-dependence  for  this  rate  and  thus  provide  a  chal¬ 
lenge  to  continuing  theoretical  efforts  to  understand  this  problem. 

More  importantly,  however,  the  combination  of  the  recent  low  tempera¬ 
ture  results  and  the  current  measurements  allows  confident  specifica¬ 
tion  of  the  Na-  0  rate  at  the  temperatures  of  the  upper  atmosphere. 
Recent  studies  show  that  with  the  use  of  these  new  rates,19-20  the 
calculated  vibrational  temperature  of  Na  in  the  E  region  is  much  lower 
than  was  previoulsy  predicted21  and  is  in  agreement  with  experiment. 


The  measured  value  of  pt^  q,  is  approximately  0.04  atm-p,s  over 

the  temperature  range  from  1200  to  3000  K.  This  is  approximately  the 

same  as  earlier  measurements  of  pt_  _  aa  and  is  approximately  25  times 

ua-u, 

more  rapid  than  the  Na-0  rate  in  the  same  temperature  range.  This 
greater  rapidity  could  result  from  either  stronger  chemical  effects  in 
the  CO-O  and  0a-0  processes  or  the  possibility  of  atom  exchange  in 
these  cases. 
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C.  Conclusions 

At  the  start  of  this  research  program,  N2  vibrational  temperatures 

21 

as  high  as  3000  K  were  predicted  for  the  E  region.  This  high  vibra¬ 
tional  temperature  resulted  from  the  balance  between  efficient  produc¬ 
tion  by  the  reactions 

N(4S)  +  NO  -*  N2(v)  +0(3P)  (2) 

0(1D)  +  Na  -  N2(v)  +  0(3P)  (3) 

and  comparatively  inefficient  removal  by  energy  transfer  to  C02; 
that  is, 

Na(v=l)  +  C02(00°0)  -*  N2(v=0)  +  C02(00°l)  .  (1) 

During  the  course  of  this  contract,  McNeal  showed  that  V-T  relaxa¬ 
tion  of  Na(v=l)  by  0(3P)  was  efficient  at  300  K  and  should  be  considered 
in  determining  the  vibrational  temperature  of  N2  in  the  E  region.1 
This  was  in  direct  conflict  with  the  theoretical  predictions  of  Fisher 
and  Bauer.6  The  measurements  performed  under  the  current  contract  ex¬ 
tend  over  the  temperature  range  from  1200  to  3000  K,  and  confirm  that  re¬ 
laxation  of  N2(v=l)  by  0(3P)  is  indeed  efficient.  We  have  also  estab¬ 
lished  that  process  (2)  is  an  efficient  source  of  vibrationally  excited 
nitrogen.  The  net  effect,  however,  is  to  expect  large  reduction  in 
the  nitrogen  vibrational  temperature  in  the  E  region  (to  ~  1000  K)  and 
to  greatly  reduce  the  sensitivity  of  the  calculated  vibrational  tem¬ 
perature  to  the  yield  of  vibrational  energy  from  the  0(JD)  +  N2  reac¬ 
tion. 


12 


REFERENCES 


1.  R.  J.  McNeal,  H.  E.  Whitson,  Jr.,  and  G.  R.  Cook,  Chem.  Phys.  Letters 
16,  507  (1972);  Trans.  Aner.  Geophys.  Union  (to  be  published). 

2.  W.  D.  Breshears  and  P.  F.  Bird,  J.  Chem.  Phys.  48,  4768  (1968). 

3.  J.  E.  Morgan,  L.  F.  Phillips,  and  H.  I.  Schiff,  Disc.  Faraday  Soc. 

33,  118  (1962). 

4.  L.  F.  Phillips  and  H.  I.  Schiff,  J.  Chem.  Phys.  36,  3283  (1962). 

5.  J.  E.  Morgan  and  H.  I.  Schiff,  Can.  J.  Chem.  41,  903  (1963). 

6.  E.  R.  Fisher  and  E.  Bauer,  J.  Chem.  Phys.  57,  1966  (1972). 

7.  G.  Black  and  D.  J.  Eckstrom,  "Excitation  and  Deexcitation  of  Vibra¬ 
tion  in  Na  by  Oxygen  Atoms, "  Final  Report,  Contract  DAHC04-70-C-0036, 
Stanford  Research  Institute,  Menlo  Park,  Calif.  (October  1971). 

8.  J.  T.  Herron,  J.  Chem.  Phys.  35,  1138,  1961;  J.  Res.  Nat.  Bur. 

Stand.  65A,  411  (1961). 

9.  L.  F.  Phillips  and  H.  I.  Schiff,  J.  Chem.  Phys.  36,  1509  (1962). 

10.  J.  D.  Lambert,  "Atomic  and  Molecular  Processes,"  ed.  by  D.  R.  Bates, 
p.  783,  1962. 

11.  A.  Dalgarno,  Planet.  Space  Sci.  10,  19  (1963). 

12.  C.  C.  Rankin  and  J.  C.  Light,  J.  Chem.  Phys.  46,  1305  (1967). 

13.  A.  L.  Schmeltekopf,  F.  C.  Fehsenfeld,  G.  I.  Sulman,  and  E.  E.  Ferguson, 

Planet.  Space  Sci.  15,  401  (1967). 

14.  D.  Rapp  and  P.  Englander-Golden,  J.  Chem.  Phys.  40,  573  (1964);  40, 

3120  (1964). 

15.  D.  Rapp,  J.  Chem.  Phys.  43,  316  (1965). 

16.  R.  S.  Mulliken,  Phys.  Rev.  36,  611  (1930). 

17.  G.  Herzberg,  "Molecular  Spectra  and  Molecular  Structure.  I.  Spectra 
of  Diatomic  Molecules,"  p.  134  (D.  Van  Norstrand  Company,  Inc.,  1950.) 

18.  R.  S.  Hickman  and  L.  H.  Liang,  Rev.  Sci.  Inst.  43,  796  (1972). 

19.  E.  L.  Breig,  M.  E.  Brennan,  and  R.  J.  McNeal,  J.  Geophys.  Res  (1973), 
to  be  published. 


13 


\ 

V 

\ 

V 

v 

» 

20.  E.  Jamshidi,  E.  R.  Fisher,  and  R.  H.  Kmamler,  J.  Geophys.  Res.  (1973), 
to  be  published. 

21.  J.  C.  G.  Walker,  Planet  Space  Sci.  16,  321,(1968);  J.  C.  G.  Walker, 

R.  S.  Stolarski,  and  A.  F.  Nagy,  Ann.  Geophys.  25,  831,  1969. 

22.  J.  H.  Kiefer  and  R.  W.  Lutz,  Eleventh  Symp.  (int..)  on  Combustion, 

The  Combustion  Institute,  Pittsburgh,  1967. 


14 


APPENDIX  A 


MEASUREMENTS  OF  VIBRATIONALLY  EXCITED  MOLECULES 
BY  RAMAN  SCATTERING.  I.  THE  YIELD 
OF  VIBRATIONALLY  EXCITED  NITROGEN  IN  THE  REACTION 
N  +  NO  -  Na  +  0 


15 


MEASUREMENTS  OF  VIBRATIONALLY  EXCITED  MOLECULES  BY  RAMAN  SCATTERING.  I. 
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Graham  Black,  Robert  L.  Sharpless,  and  Tom  G.  Slanger 
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The  reaction  of  nitrogen  atoms  with  nitric  oxide  leads  to  vibra- 
tionally  excited  nitrogen.  In  the  presence  of  molecular  nitrogen,  V-V 
exchange  processes  rapidly  establish  a  modified  distribution  for  which 
a  measurement  of  N  (v=l)  is  adequate  to  determine  the  total  vibrational 
energy  yield  of  the  reaction.  The  measurements  of  N (v**l)  have  been 

a 

made  by  Raman  spectroscopy  using  the  intensity  of  the  Q  branch  of  the 
anti-Stokes  line  at  4382  A  (using  4880  A  Ar+  laser  radiation).  The 
system  sensitivity  is  determined  using  heated  nitrogen.  The  measure¬ 
ments  have  established  that  25  ±  3#  of  the  available  energy  (3.27  eV) 
appears  as  vibrational  energy  and  that  the  vibrationally  excited  nitro¬ 
gen  is  an  initial  product  and  does  not  result  from  collisions  of  trans- 
lationally  hot  oxygen  atoms  with  nitrogen. 


I,  INTRODUCTION 

Vibrationally  excited  nitrogen  plays  a  number  of  important  roles 
in  determining  the  nature  of  the  lower  thermosphere— its  influence  on 
the  ionospheric  composition, 1-3  the  electron  temperature,4"®  and  the 
energy  balance7  are  of  special  importance. 

In  the  undistrubed  E  region  the  two  most  important  processes  that 
produce  vibrationally  excited  nitrogen  are7 

0(*D)  +  N2  -N2(v)  +  0(3P)  (1) 

and 

* 

N(4S)  +  NO  -  N  (v)  ’+  0(3P)  .  \  (2) 

4u 

A  calculation  by  Fisher  and  Bauer®  indicated  that  less  than  5$  of  the 
available  energy  in  reaction  (l)  (1.96  eV)  is  expected  to  channel  into 
the  vibrational  levels  of  Ng.  An  experiment  to  determine  the  vibra¬ 
tional  yield  of  reaction  (l)  was  recently  carried  out  in  this  laboratory,® 
but  the  results  were  inconclusive. 10  An  efficiency  less  than  5$  makes 
reaction  (l)  a  minor  source  of  N  vibrational  energy  in  the  E  region. 

Measurements  of  reaction  (2)  a  decade  ago  by  Schiff  and  coworkers11"13 
indicated  that  about  30$  of  the  exothermicity  of  the  reaction  appears  as 
vibrational  energy.  Wray14  has  pointed  out  an  alternative  interpretation 
of  the  earlier  measurements,  which  allows  the  exothermicity  of  the  reac¬ 
tion  to  appear  as  translational  energy  of  the  product  species — prefer¬ 
entially  in  the  oxygen  atom  because  of  its  lighter  mass.  The  transla- 
tionally  hot  oxygen  atoms,  under  the  experimental  conditions,  would 
undergo  collisions  with  nitrogen.  These  collisions  would  have  a  high 

A-2 
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probability  si  exciting  vibrational  energy  in  the  nitrogen. 15  The 
validity  of  the  earlier  measurements  can  also  be  questioned  on  the 
grounds  that  the  calorimetric  probe  used  may  have  responded  to  energetic 
species  other  than  N^(v). 

The  experiments  to  be  described  below  havf  determined  the  fraction 
of  the  exothermicity  of  reaction  (2)  which  appears  as  vibrational  en¬ 
ergy,  using  Raman  scattering  as  a  specific  detector  of  vibrationally 
excited  nitrogen  molecules.  The  experimental  conditions  have  also  been 
varied  to  show  that  the  vibrationally  excited  nitrogen  is  a  primary 
product  of  reaction  (2)  and  that  any  production  of  N  (v=l)  by  collisions 
of  translationally  hot  oxygen  atoms  with  N  ,  if  it  occurs  at  all,  must 
be  of  minor  importance. 


II.  EXPERIMENTAL 


Reaction  (2)  is  fast  («  3x10" 11  cm3  molec"1  sec"1  at  room  tempera¬ 
ture16)  and  is  sufficiently  exothermic  (3.27  ev)  to  produce  N  (v)  up  to 

2 

v=12.  Whatever  the  intlal  vibrational  distribution,  it  will  be  rapidly 
degraded  by  V-V  exchange  collisions  of  the  type 

N2(v=n)  +N2(v=0)  -  N2(v<=n-n')  +  N2(v=n')  (3) 

resulting  in  an  Increased  population  of  N  (v=l).  (At  N  pressures  of 

2  2 

several  torr  this  will  require  only  a  few  milliseconds.  Loss  of  vibra¬ 
tional  energy  by  diffusion  and  wall  deactivation  is  small  in  this  time, 

and  is  taken  into  account.)  Hence  a  measurement  of  thermalized  N  (v=l) 

2 

is  sufficient  to  determine  the  fraction  of  the  exothermicity  of  reaction 


(2)  that  appears  as  vibrational  energy. 


The  principle  of  the  experiment  is  to  make  N  (vib)  by  reaction  (2) 

2 

followed  by  reaction  (3),  and  to  monitor  the  N  (v=l)  concentration  by 
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o  + 

observing  the  Raman  scattering  of  4880  A  Ar  laser  radiation.  The 
intensity  of  the  Q  branch  of  the  anti-Stokes  scattering  (v=l  -♦  v=0)  at 

O 

4382  A  was  measured  because  the  anti-Stokes  spectra  are  free  of  Raman 
scattering  from  the  Ng(v=o)  level  (whereas  the  Stokes  transitions  from 
v=l  -*  v=2  lie  in  the  rotational  structure  of  the  Stokes  lines  from 
v=0  -*  v=l) . 

The  experimental  arrangement  is  shown  schematically  in  Fig.  1. 
Molecular  nitrogen  is  partially  dissociated  by  a  microwave  discharge  in 
a  flow  system  having  a  linear  velocity  of  400  cm  sec-1.  The  nitrogen 
atoms  can  be  titrated  with  NO  at  any  of  six  fixed  inlet  jets  at  8,  15, 
30,  42,  65,  and  100  cm  from  the  pyrex  cell  (referred  to  as  inlets  1 
through  6,  respectively).  The  output  of  an  Ar+  ion  laser  (Coherent 
Radiation  Model  52B)  traverses  the  cell  and  is  focused  to  a  point 
within  it.  This  focal  point  is  imaged  on  the  slit  of  a  spectrometer 
(Spex  Model  1402)  by  a  3econd  optical  system.  The  light  passing  through 
the  exit  slit  of  this  spectrometer  is  imaged  on  the  photocathode  of  a 
6256S  photomultiplier.  Its  output  is  measured  in  a  pulse-counting  sys¬ 
tem  and  then  printed  out  after  a  predetermined  counting  time  (governed 
by  the  setting  of  tne  controller). 

Initial  experiments  were  made  to  determine  the  magnitude  of  the 

interference  from  emissions  of  the  titrated  afterglow.  The  nature  of 

this  problem  can  be  seen  in  Fig.  2,  which  shows  the  light  signal  at 
0 

4382  A  from  the  afterglow  as  a  function  of  NO  addition.  The  signal 
reaches  a  minimum  of  ~  10  counts  sec"1  (which  includes  a  dark  count  of 
~5  sec”1)  when  the  NO  is  just  sufficient  to  remove  the  N  atoms  by  reac¬ 
tion  (2)  (called  "null").  Further  addition  of  NO  generates  the  NO 
continuum  by  the  reaction  of  the  excess  NO  with  the  oxygen  atoms,  that 


followed  by 


NO  -NO  +  hv  (3875-14,000  A)  .  (5) 

A  2 

Thus  when  NO  is  in  excess,  radiation  from  (5)  is  transmitted  by  the 
spectrometer  at  4382  A. 

When  insufficient  NO  is  added,  the  signal  arises  from  the  genera¬ 
tion  of  NO  3  bands  by  the  chemiluminescent  recombination  of  oxygen  and 
nitrogen  atoms,  that  is, 

N  +  0  +  M  -  NO  (Ban  )  +  M  (6) 

r 

N0(Banr)  -4  NO  (xanr)  +  hv(g  bands)  .  (7) 

The  light  being  transmitted  by  the  spectrometer  at  4382  A  is  in  the 
tail  of  the  (0,13)  band,  whose  head  is  at  4309.7  A  (the  bands  are  de¬ 
graded  to  longer  wavelengths).  (The  bandpass  of  the  spectrometer 
(~  1*1  A)  was  such  as  not  to  transmit  any  of  the  short  wavelength  head 
of  the  (1,14)  band  at  4385*7  A.) 

The  trace  on  Fig.  2  clearly  demonstrates  the  need  to  maintain  the 
NO  addition  as  close  to  the  null  point  as  possible.  To  do  this  the 
valve  controlling  the  NO  addition  was  coupled  to  a  servo-controlled 
motor.  The  inputs  to  the  motor  controller  were  derived  from  two  1P21 
photomultipliers  viewing  the  afterglow  downstream  of  the  pyrex  cell. 

One  photomultiplier  "sensed"  the  continuum  through  a  5622  A  interference 
filter  (10  X  half-width),  while  the  other  viewed  the  afterglow  through 
a  Kodak  18A  filter  and  was  therefore  sensitive  only  to  the  NO  g  bands. 
With  this  arrangement,  it  was  possible  to  maintain  the  NO  addition 
within  1$  of  null.  Therefore,  for  Raman  measurements  on  the  titrated 
afterglow,  the  servo-controlled  motor  was  always  used  to  regulate  the 
NO  addition. 
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Since  reaction  (2)  is  the  well-know  titration  reaction  for  nitrogen 
atoms  the  partial  pressure  of  nitrogen  atoms  being  reacted  with  MO  was 
determined  in  the  following  manner.  With  the  microwave  discharge  on  the 
servo-controlled  motor  was  used  to  regulate  the  NO  addition  to  the  null 
position  at  which  the  partial  pressure  of  NO  being  added  equals  the  par¬ 
tial  pressure  of  nitrogen  atoms  in  the  afterglow.  The  NO  partial  pressure 
being  added  was  then  determined  by  turning  both  the  servo-controlled  motor 
and  the  microwave  discharge  off  and  determining  the  change  in  the  system 
pressure  when  the  NO  was  turned  off,  using  a  differential  pressure  dia¬ 
phragm  transducer  attached  to  the  vacuum  system  downstream  of  the  Pyrex 
cell.  This  pressure  change  was  too  small  for  accurate  measurement  with 
100$  NO  so  an  analyzed  mixture  of  10$  NO  in  helium  was  used  for  the  mea¬ 
surements. 

The  measurements  were  made  with  an  EMI  6256S  phototube  with  20-dB 
attenuation  inserted  in  the  line  between  the  photomultiplier  base  and 
the  photon  counter,  to  provide  the  greatest  freedom  from  rf  pickup. 

The  best  signal-to-dark  count  was  obtained  with  1200  V  on  the  photo¬ 
multiplier,  which  resulted  in  a  counting  efficiency  75$  of  the  plateau 
maximum  and  a  dark  count  of  ~  5  sec**1. 

Before  making  any  Raman  measurements,  the  positions  of  lenses 
and  L  were  adjusted  to  provide  the  biggest  Rayleigh-scattered  light 
signal  (with  N  at  600  torr  in  the  cell)  using  the  4579  A  laser  line 
(the  shortest  wavelength  laser  line  and  hence  that  closest  in  wavelength 
to  the  position  of  the  anti-Stokes  line  at  4382  A).  This  focus  was 
periodically  checked  but  was  found  to  require  minimum  adjustment.  The 
wavelength  error  in  the  readout  of  the  spectrometer  was  periodically 
checked  with  the  4358.3  A  line  of  Hg.  It  did  not  change  during  the 
course  of  these  measurements. 
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In  addition  to  several  photomultiplier  counts  per  second  originat¬ 
ing  from  the  titrated  afterglow,  a  small  signal  (<  0.03  counts/sec)  at 
4382  A  was  also  produced  by  the  laser  (arising  from  the  4379.7  A  argon 
line  emitted  by  the  discharge  in  the  laser  cavity).  Because  of  these 
two  sources  of  background,  the  Raman  signal  had  to  be  extracted  from 
the  counts  obtained  in  four  consecutive  counting  periods  of  equal 
length.  From  such  a  sequence,  which  we  define  as  an  "experiment,”  the 
Raman  signal  is  obtained  from 

Raman  counts  =  -  Cg  ,  (8) 

where  is  the  count  obtained  with  the  microwave  discharge  and  laser 

on  (and  hence  includes  the  Raman  signal  together  with  the  scattered 

light  from  both  sources  and  the  dark  count):  C„  is  the  photomultiplier 

4 

dark  count  (microwave  discharge  and  laser  off);  is  the  count  with 

only  the  microwave  discharge  on;  and  C  is  the  count  with  only  the  laser 

3 

on.  For  most  of  the  experiments,  the  counting  period  was  six  minutes. 

It  was  soon  discovered  that,  even  without  titrating  the  N  atoms,  the 
N  (v»l)  present  in  the  afterglow  was  sufficient  to  produce  a  large  Raman 
signal  (increasing  from  O.S  to  20  counts/sec  as  the  N  pressure  was  raised 
from  1.5  to  12  torr).  To  increase  the  detectability  of  the  N  (v=l)  ex- 

4d 

pected  from  the  N  +  NO  reaction,  it  was  decided  to  remove  or  at  least 
minimize  the  N  (v«*l)  originating  in  the  discharge.  This  was  done  by  in- 
sorting  a  glass-wool  plug  between  the  discharge  and  the  NO-inlet  positions. 
The  logic  behind  this  approach,  which  has  been  used  before11'19  is  that 
the  deactivation  of  N  (v«l)  on  pyrex  surfaces  (v  =  4.5  x  10“4)  is  more 

d 

efficient  than  the  removal  of  nitrogen  atoms  (Y  a  2  x  10-8). 17-18  The 
first  two  trials  of  a  glass-wool  plug  confirmed  this  situation,  reducing 
the  N _(v*l)  in  N  at  6  torr  by  factors  of  3  and  5,  while  reducing  the 
nitrogen  atoms  by  10  and  25$,  respectively.  Several  large  glass-wool 

zz 
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plugs  were  tried— the  largest  decreasing  the  N  (v=l)  in  K  at  6  torr  by 

2  2 

a  factor  of  30  (to  ~  0.3  counts/sec)  while  reducing  the  nitrogen  atoms  by 

two  (to  ~8  mtorr).  Host  of  the  measurements  to  be  described  were  made 

under  these  conditions.  The  final  measurements  were  made  with  a  5$  N 

in  argon  mixture  at  10  torr  (with  the  largest  pyrex  wool  plug  in  place ), 

which  reduced  the  N (v=l)  signal  with  no  titration  by  another  factor  of 
2 

five  (to  ~  0.06  counts/sec)  while  only  reducing  the  N  atoms  by  another 
40#  (to  ~  5  mtorr). 

Because  it  took  many  minutes  for  the  microwave  discharge  to  stabilize, 
it  was  found  preferable  to  leave  the  microwave  discharge  on  all  the  time 
and  to  evaluate  the  Raman  counts  plus  a  small  signal  due  to  the  laser 
scattered  light  from: 

Raman  counts  +  laser  scattered  light  =  «  S  *  (9) 


This  was  determined  for  both  the  untitrated  afterglow  (=»  S  ^  J 

untitrated 

and  the  titrated  afterglow  when  the  titration  was  performed  at  inlets 
no.  1,  2,  3,  4,  and  sometimes  5.  The  titrated  signals  were  then  plotted 
as  a  function  of  distance  from  the  titration  point  to  the  cell  and  extra¬ 
polated  to  zero  distance  to  determine  the  signal  S  ,  tnat  would  have  been 

o 

produced  if  there  had  been  no  loss  of  N  (vel)  between  the  titration  point 

2 

and  the  cell.  From  this  signal,  S  .. .  .  .  was  subtracted  to  obtain  the 

untitrated 

signal  produced  by  reaction  (2),  (s  ).  The  subtraction  removes  the 

reaction 

small  signal  from  the  scattered  laser  light  which  is  present  in  both  S 

o 

and  S  • 

untitrated 
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III.  CALIBRATION 


To  convert  the  net  measured  Raman  signal  into  a  partial  pressure  of 
N  (v«l)  at  297*$  we  measured  the  intensity  of  the  anti-Stokes  line  at 

A 

o 

4382  A  with  nitrogen  at  atmospheric  pressure  in  the  cell.  Table  I  shows 
the  data  obtained  using  a  chromel-alumel  thermocouple  to  measure  the 
nitrogen  temperature.  The  last  column  corrects  the  actual  pressure  of 
N  (v=l)  to  the  equivalent  pressure  at  297°K.  This  is  shown  plotted 

A 

against  photomultiplier  counting  rate  in  Fig.  3.  The  horizontal  error 
bars  show  the  effect  of  an  uncertainty  of  ±  2*K  in  the  gas  temperature. 

The  straight  line  drawn  gives  a  sensitivity  of  41  counts/sec  [torr  N  (v=l)]~ 
This  experiment  was  repeated  several  times  during  the  measurement  pro¬ 
gram  and  was  always  found  to  lie  within  a  range  41  i  3  counts/sec 
^torr  N  (v=l)]-1. 

A 

Because  thermal  conduction  might  affect  the  thermocouple  readings, 
we  also  determined  the  gas  temperature  from  the  intensity  variation  of 
the  rotational  Raman  spectrum.  This  rotational  Raman  spectrum  is  of  high 
quality,  and  provides  a  reliable  temperature  measurement.  The  intensity 
formula  is90 


I  =  f(v)  v3s  exp(-BJ(j+l)hc/kT)  ,  (10) 

U 

where  f(v)  incorporates  the  possible  frequency  dependence  of  the  counting 
efficiency,  S  is  the  line  strength,  and  the  v3  variation  must  not  be 
ignored  if  the  correct  T  is  to  be  obtained.  The  Stokes  lines  progress 
to  larger  wavelengths  with  J,  and  the  anti-Stokei<  lines  progress  to 
smaller  wavelengths,  but  they  give  a  common  temperature  to  within  i  1  K 
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when  f ( v)  is  assumed  constant.  This  indicates  that  the  frequency  depen¬ 
dence  of  the  counting  efficiency  is  negligible  over  the  range  (~  100  A) 
of  the  rotational  spectrum.  The  resulting  rotational  temperature  agrees 
with  the  thermocouple  reading  at  room  temperature,  and  the  thermocouple 

O  0 

reading  was  about  4  K  low  when  the  rotational  temperature  was  377  K.  We 
have  therefore  made  a  linear  correction  to  the  thermocouple  reading, 
varying  from  0 <f>  at  297  °K  to  +  1#  at  373 °K.  These  corrected  temperatures 
are  used  in  Table  I. 


TABLE  1.  Photomultiplier  Count  Rate  and  Pressure 
of  N  (v«l)  in  N  at  Various  Temperatures 

6  A 


Photomultiplier 
Temperature  Signal 

(T  k)  (counts/sec) 

[n2(v*D] 

[n2(v=o)] 

[Nz(v=l)] 

(Tor r  in  N 

A 

at  760  Torr) 

[N2(v«1)]x 

297  m 
-vr-  Torr 

T  K 

297 

5.44 

1.25  x  10-6 

0.0095 

0.0095 

323 

5.89 

3.10  x  10“5 

0.024 

0.022 

350 

6.88 

6.90  X  10"8 

0.053 

0.045 

364 

7.49 

9.98  x  10"8 

0.076 

0.062 

379 

8.51 

1.44  x  10-4 

0. 109 

0.086 
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IV.  RESULTS 


The  Raman  counts  obtained  from  the  titration  were  about  1  sec-1,  and 
this  signal  had  to  be  measured  in  the  presence  of  a  dark  count  of  5  sec-1 
and  several  counts  per  second  produced  by  light  from  the  titrated  after¬ 
glow.  Therefore  a  large  number  of  experiments  had  to  be  performed  at 
each  titration  point  to  determine  the  signal  with  sufficient  accuracy. 

Table  II  is  an  example  of  the  data  obtained  in  approximately  six  hours  at 
one  titration  point.  This  data  point  and  the  ones  determined  when 
titrating  at  inlets  no.  1,  3,  4  and  5  are  shown  in  Fig.  4.  It  can  be  seen 
that,  with  increasing  distance,  the  signals  obtained  for  the  titrated 
afterglow  appear  to  be  asymptotic  to  the  value  obtained  for  the  untitrated 
afterglow.  This  is  because  the  microwave  discharge  produces  so  much  more 

N  (v=l)  than  the  titration  that  as  the  titration  point  moves  further  from 
2 

the  Pyrex  cell  and  closer  to  the  discharge,  the  N2(v=l)  produced  by  the 
titration  reaction  becomes  a  smaller  and  smaller  fraction  of  the  N2(v=l) 
surviving  from  the  discharge.  The  error  bars  are  the  estimated  error  of 
the  mean,  p-,  determined  from 

*  v 0 


a 


where  a  is  the  standard  deviation  of  the  individual  x  values  and  N  is 
x 

the  number  os  such  values.  At  various  times  during  the  measurement  the 
NO  was  turned  off  and  the  Raman  signal  from  the  untitrated  afterglow  was 
determined,  to  ensure  that  neither  the  optical  focus  nor  the  aiterglow 
condition  had  changed. 

The  result  obtained  in  four  separate  series  of  measvrements  are 
shown  in  Table  III.  They  have  been  arranged  in  order  of  decreasing  Raman 
signals  from  the  untitrated  afterglow.  The  Raman  signal  produced  by  the 


r 


▼ 


TABLE  II.  Signal  Obtained  by  Titrating  8  mtorr 
of  N  Atoms  at  Inlet  No.  2 

(Ng  =  6  torr) 


1 


Expt. 

No. 

Laser  off  Laser  on 

C,  Counts  C„  Counts 
1  2 

C2-°l 

Counts 

Expt. 

No. 

Laser  off 
C^  Counts 

Laser  on 

C  Counts 
2 

Vci 

Counts 

3884 

4196 

312 

15 

3958 

4189 

231 

3875 

4329 

454 

16 

4060 

4235 

175 

4028 

4197 

169 

17 

3939 

4250 

311 

3694 

4204 

510 

18 

3890 

4384 

494 

3949 

4116 

167 

19 

3959 

4213 

254 

3948 

4106 

15'J 

20 

4053 

4248 

195 

3848 

4044 

196 

21 

4003 

4219 

216 

8 

3921 

4095 

174 

22 

3903 

4381 

478 

9 

3970 

4393 

423 

23 

4029 

4370 

341 

10 

3851 

4191 

340 

24 

4063 

4321 

258 

11 

4000 

4202 

202 

25 

3949 

4262 

313 

12 

3965 

4333 

368 

26 

3998 

4304 

306 

13 

3784 

4215 

431 

27 

4047 

4349 

302 

14 

393? 

4289 

357 

28 

4012 

4526 

514 

Note:  Counts  in  a  6  minute  period.  Signal 
from  untitrated  afterglow  =  10?  ±  11 


Mean  x  =  309 

Number  of  measurements,  N  =  28 

a  =  115 
x 

a 

x 

Estimated  error  of  mean  p-  =  — —  =  22 

*  >/* 
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TABLE  III.  Efficiency  with  which  the  Exothermicity  of  the 
Reaction  of  Nitrogen  Atoms  with  Nitric  Oxide  is 
Converted  into  Vibrational  Energy. 


Gas 

Pressure 

S 

untitrated 

(counts/sec 

S 

0 

)  (counts /sec) 

s  ■*  s 

0  untitrated 

=  Reaction 
(counts/sec) 

N2(v=1)  N 

(mtorr)(mtorr) 

F 

N2 

6  Torr 

1.17 

2.18 

1.01 

25 

11 

0.20 

N2 

6  Torr 

0.30 

1.11 

0.81 

20 

8.0 

0.22  | 

"2 

6  Torr 

0.29 

1.12 

0.83 

20 

7.5 

0.24 

5*  V 

argon 

10  Torr 

0.056 

0.68 

0.62 

15 

5.5 

0.25 

reaction  is  shown  in  column  5  and  is  then  converted  to  a  N  ( v=l)  pressure 

A 

using  the  sensitivity  of  41  counts/sec  [torr  N  (vsl)]"1  —  this  is  shown 
in  column  6.  Column  7  shows  the  N  atom  pressure.  The  efficiency  with 
which  the  exothermicity  of  the  reaction  is  converted  into  vibrational 
energy  is  shown  in  the  final  column.  This  efficiency,  F,  is  calculated 
from 


Vv=i) 
11,3  p» 


(12) 
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where  ^  ^  is  the  partial  pressure  of  N0(v=l)  produced  by  reacting 

a  partial  pressure,  P  ,  of  nitrogen  atoms  with  NO  and  the  factor  11.3 

is  introduced  because  the  exothermicity  of  the  reaction  (3.27  eV)  is 

11.3  times  the  vibrational  energy  of  N  (v=l)  [0.289  eV]. 

2 


The  value  of  F  determined  with  5$  N^  in  argon  at  10  torr  has  been 

increased  by  6$  above  that  calculated  from  Eq.  (12)  to  allow  for  the 

Boltzmann  population  in  the  v=2  level  ( for  a  vibrational  temperature  of 
0  . 

960  K).  This  is  a  negligible  correction  in  the  other  three  sets  of  ex¬ 


periments  with  N  at  6  torr,  because  in  these  cases  the  N  vibrational 
2  2 

temperature  was  much  smaller  610  K). 


V.  DISCUSSION 

The  value  of  F  in  Table  III  tends  to  increase  as  the  background 

(S  )  is  reduced.  It  was  originally  thought  that  this  might  re¬ 

untitrated 

fleet  an  increased  rate  of  decay  of  N  (v«*l)  in  the  presence  of  oxygen 

atoms  (produced  by  reaction  (2))  compared  to  when  only  nitrogen  atoms 

are  present.  Part  of  this  increased  rate  of  decay  might  involve  the 

V-T  relaxation  of  N  (v==l)  by  0(3P),  which  has  recently21  been  shown  to 
2 

have  a  rate  of  3.5  x  10"15  cm3  molec-1  sec-1  and  part  might  arise  from 

an  increased  catalytic  efficiency  of  the  walls.  However,  extrapolating 

the  measurements  to  zero  distance  should  remove  any  homogeneous  or  wall 

decay  between  the  inlet  jet  and  the  cell.  The  possibility  remains  that 

the  walls  of  the  inlet  jet  itself  or  hot  0(3P)  produced  in  the  vicinity 

of  the  inlet  jet  by  reaction  (2)  could  deactivate  N (v=l).  These  effects 

2 

would  not  be  allowed  for  by  the  extrapolation  to  zero  distance.  It  is 
for  these  reasons  that  the  values  obtained  with  the  lowest  background 
from  the  discharge  are  to  be  most  trusted.  With  this  in  mind  and  allow¬ 
ing  for  possible  errors  in  the  extrapolation  to  determine  and  in  mea¬ 
suring  the  N  atom  concentration  and  the  sensitivity  to  N  (v=l),  we 

2 
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recommend  a  value  of  0.25  ±  0.03  for  the  efficiency  with  which  the  exo- 
thermicity  of  reaction  (2)  appear  as  N  vibrational  energy. 

It  should  also  be  noted  that  the  last  value  of  F  in  Table  III  was 
obtained  with  only  0.5  torr  of  N  (+9.5  torr  of  argon)  and  is  in  agree- 

<u 

ment  with  the  measurements  at  N  =  6  torr.  This  requires  that  0.5  torr 

2 

of  N  be  more  than  sufficient  to  Boltzmannize  the  initial  distribution 
2 

of  vibrational  levels  in  a  flow  time  of  0.02  sec  (8  cm  at  400  cm  sec"1). 

This  corresponds  to  a  rate  coefficient  >  3  x  10“18  cm3  molec"1  sec-1 
which  is  much  faster  than  the  value  of  3.5  x  10"1G  cm3  molec-1  sec-1 
found  by  Phillips  and  Schiff12  for  colllsional  deactivation  of  N  (v)  to 
vibrational  levels  below  the  fourth.  The  reason  for  this  discrepancy 
is  not  apparent  although  Wray14  has  questioned  the  low  rate  found  in  the 
earlier  work.  [A  possible  explanation,  which  seems  very  unlikely  but 
cannot  be  excluded,  is  that  decomposition  of  ozone  is  possible  by  step¬ 
wise  excitation  of  ozone  by  levels  below  v  =  4.]  » 

The  measurement  with  5$  N  in  argon  shows  that  the  vibrationally 
excited  nitrogen  is  a  product  of  the  reaction  (2)  and  is  not  subse¬ 
quently  excited  by  collisions  of  translationally  hot  oxygen  atoms, 
produced  by  reaction  (2),  with  N  .  If  this  latter  process  were  involved 
the  large  excess  of  argon'  over  nitrogen  in  the  latter  measurement  would 
thermalize  the  hot  oxygen  atoms  and  prevent  vibrational  excitation  of  Ng. 
Hence  Wray's  criticism  of  the  earlier  work11-13  on  this  reaction  does 
not  appear  to  be  valid.  In  fact  the  value  obtained  by  Schiff  and  his 
coworkers  (F  =  0.27  ±  0.06)  is  in  agreement  with  our  value  despite  the 
nonspecific  nature  of  their  calorimetric  detector. 
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FIGURE  CAPTIONS 


Fig.  1  Experimental  arrangement  for  Raman  spectroscopy.  The  lens 

L  ,  which  has  an  aperture  ratio  of  f/0.95  and  a  focal  length 
of  5  cm,  collects  the  Raman-scattered  light  over  a  large 
solid  angle.  The  titrated  afterglow  enters  the  pyrex  cell 
immediately  above  the  focus  of  the  laser  beam. 

O 

Fig.  2  Light  signal  at  4382  A  from  microwave-discharged  N^  at 
3  torr  as  a  function  of  NO  addition  close  to  the  equiv- 

Q 

alence  point  ("null").  Spectrometer  bandpass  ~1*1  A. 
Titration  performed  100  cm  (~  0*25  second  flow  time) 
from  cell. 

Fig.  3  Determination  of  the  experimental  sensitivity  to  N  (v=l) 

a 

using  heated  N  . 

2 

Fig.  4  Raman  signal  produced  when  8  mtorr  of  nitre* an  atoms  are 
reacted  with  NO  at  various  distances  from  the  center  of 
the  pyrex  cell.  Ng  =  6  torr.  With  no  titration  a  Raman 
signal  of  0*30  counts /sec  is  present  in  the  afterglow. 
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ABSTRACT 


We  have  studied  the  effects  of  oxygen  atoms  on  the  vibrational 
relaxation  of  shock -heated  Na  over  the  temperature  range  from  1200  to 
3000  K,  using  the  00  tracer  technique.  The  measured  relaxation  times 
PT  _  vary  from  0.8  to  2  atm-p,s  in  this  range  with  an  uncertainty  of 

Nj|"U 

±  50$,  which  is  in  excellent  agreement  with  existing  high-temperature 
shock  tube  and  low- temperature  flow  tube  results.  The  weak  temperature 
dependence  of  all  of  these  results  is  in  marked  contrast  to  theoretical 
predictions.  However,  the  relatively  fast  relaxation  rates  at  low  tem¬ 
peratures  verified  by  these  measurements  allow  a  satisfactory  explana¬ 
tion  of  the  observed  vibrational  temperature  of  Na  in  the  upper  atmosphere. 


Some  limited  measurements  of  the  relaxation  of  CO  by  oxygen  atoms 
yield  values  of  of  approximately  0.04  atm-n,s.  This  is  approxi¬ 
mately  the  same  as  pt  ,  and  probably  is  indicative  of  the  importance 
of  atom-exchange  processes  in  both  of  these  reactions. 
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I  INTRODUCTION 


A  number  of  experiments  have  been  reported  in  which  the  vibrational 
relaxation  of  a  diatomic  gas  by  a  chemically  active  atomic  species  was 
studied.  These  include  shock  tube  studies  of  the  effects  of  0(3P)  on  03 
by  Kiefer  and  Lutz,1  of  0(3P)  on  N2  by  Breshears  and  Bird,3  of  H  on  CO 
and  of  Fe  on  CO  by  von  Rosenberg,  et  al.^4  Very  recently,  McNeal, 
Whitson,  and  Cook5  have  reported  low-temperature  measurements  in  a  flow¬ 
ing  afterglow  of  the  effect  of  0(3P)  on  Na.  In  all  cases,  the  relaxa¬ 
tion  times  for  these  diatom-atom  collisions  are  considerably  shorter 
than  would  be  predicted  from  correlation::  based  on  reduced  mass  (e.g. , 
Millikan  and  White6);  furthermore,  the  relaxation  times  generally  show 
very  little  temperature-dependence. 

For  the  0-0a  and  0-Na  cases,  and  speculatively  for  the  Fe-CO  case, 
the  rapid  vibrational  energy  transfer  has  been  attributed  to  a  "chemical" 
effect  since  in  each  case  there  can  be  an  attractive  chemical  bond. 

Fisher  and  Bauer,7  using  a  curve-crossing  model,  have  made  calculations 
for  the  0-Na  system  that  predict  the  high-temperature  results  of  Breshears 
and  Bird.  Their  assumed  potential  for  the  interaction  of  Na  with  0  gives 
an  activation  energy  of  0.8  eV;  as  a  result,  their  calculated  results  show 
a  very  strong  temperature-dependence.  Therefore,  their  predicted  room 
temperature  relaxation  rate  is  more  than  8  orders  of  magnitude  slower  than 
the  experimental  results  of  McNeal,  et  al. 

Breshears  and  Bird’s  measurements  with  the  laser  schlieren  technique 
covered  the  temperature  range  from  3000  to  4500  K.  The  lower  temperature 
limit  was  imposed  by  the  onset  of  large  random  fluctuations  in  the 
schlieren  signal  that  were  probably  caused  by  turbulence  in  the  shock 
tube  boundary  layers.  We  report  here  the  results  of  a  program  undertaken 
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to  extend  the  lower  temperature  limit  of  the  shock  tube  measurements  of 
the  0-N2  system  to  approximately  1000  K  by  the  use  of  the  CO  infrared 
tracer  technique.8  This  method  has  been  used  to  measure  the  relaxation 
of  N2  and  of  N2-02  mixtures  for  temperatures  as  low  as  800  K,8-10 
which  is  well  below  the  limit  for  interferometry  or  laser  schlieren  mea¬ 
surements  for  these  gases.  The  CO  tracer  method  has  a  shortcoming,  how¬ 
ever,  in  that  it  introduces  additional  species  into  the  gas  mixture, 
namely  CO  and  CO^  whose  effects  on  the  Na-0  system  must  be  considered. 
Also,  as  described  later,  components  of  IR  radiation  in  addition  to  the 
desired  CO  signal  can  possibly  interfere  with  the  measurements. 

As  in  the  other  shock  tube  studies  of  the  effects  of  0  atoms,  the 
ground  state  0(3P)  atoms  are  obtained  from  the  shock -heated  thermal  dis¬ 
sociation  of  03  added  to  the  test  gas  mixture.  The  primary  reactions  of 
interest  are 


03  +  M  -  0  +  02  +  M  ,  (l) 

0  +  03  -  02  +  02  (2) 

and 

Oa  +  M  s±  0  +  0  +  M  ,  (3) 

where  M  is  any  collision  partner,  predominantly  N2  .  At  high  tempera¬ 
tures,  the  dissociation  of  03  is  dominated  by  reaction  (l),  and  there 
is  virtually  a  one-to-one  conversion  of  03  to  0  atoms.  At  lower  tem¬ 
peratures,  however,  reaction  (2)  becomes  increasingly  important  so  that 
relatively  fewer  0  atoms  are  produced.  The  rates  for  these  processes, 
their  uncertainties,  and  the  fractional  production  of  0  atoms  are  dis¬ 
cussed  in  the  Appendix . 

The  decomposition  of  0,  behind  the  shock  wave  by  reactions  (l)  and 
(2)  is  accomplished  in  less  than  1  us  for  the  conditions  of  these  tests, 
while  reaction  (3)  is  much  slower  than  vibrational  relaxation  processes 
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and  is  therefore  unimportant.  Three  other  reactions  of  possible  im¬ 
portance  are: 


and 


C0+0+M^C0a+M  , 

CO  +  02  2  C02  +  0 

0  +  Na  st  NO  +  N 


(4) 

(5) 

(e) 


Again,  an  examination  of  the  rates  for  these  processes  indicates  that 
they  are  slow  compared  with  vibrational  relaxation  processes.  There¬ 
fore,  we  can  assume  that  the  concentration  of  0  atoms  is  constant  during 
the  time  interval  of  interest. 
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II  EXPERIMENTAL 


The  shock  tube  used  in  this  study  is  a  stainless  steel  tube  5  cm 
in  diameter  by  4  m  long,  operated  with  cold  Ha  pressure  drive.  The 
test  station  40  cm  from  the  endwall  was  fitted  with  a  sapphire  infrared 
window  and  with  two  quartz  windows  for  UV  absorption  measurements.  An 
InSb  photovoltaic  detector  was  focussed  on  the  center  of  the  shock  tube 
using  an  IRTRAN  II  lens.  The  infrared  optical  alignment  was  accomplished 
by  maximizing  the  chopped  radiation  signal  from  a  small  thermal  source 
placed  at  the  shock  tube  axis.  A  narrow-bandpass  filter  (5.05  to  5.15  p,) 
was  used  in  all  tests.  In  most  of  the  experiments,  the  detector  was 
used  in  the  voltage  mode,  and  the  signal  amplified  by  an  operational  am¬ 
plifier  having  a  gain  of  50.  The  response  time  of  this  detector-amplifier 
configuration  was  determined  to  be  1.6  (jus  when  an  IR  LED  and  an  HP  8002A 
pulse  generator  with  10  ns  risetime  were  used.  Late  in  the  program,  an 
automatically  biasing  current  mode  amplifier  was  obtained.  The  response 
time  was  then  reduced  to  less  than  1  For  all  the  data  discussed 
here,  the  intensity  risetime  was  several  times  the  detector  response 
time. 

The  Na-C0  and  Na-C0-0a  gases  used  in  the  study  were  premixed  by 
Matheson  and  were  of  at  least  99.995$  purity.  Ozone  was  made  by  dis¬ 
charging  0a  and  condensing  the  resulting  0,  at  liquid  Na  temperatures. 

The  Na-C0-03  mixtures  were  then  made  by  using  the  premixed  Na-C0  after 
passage  through  a  Millikan-type  cold  trap, 1 1  with  pressures  measured  on 
calibrated  Wallace  and  Tiernan  gages.  These  test  samples  were  stored 
at  dry  ice  temperatures  to  reduce  decomposition  of  the  0g;  samples  were 
generally  made  daily.  The  actual  concentration  of  0,  in  each  test  gas 

O 

fill  was  measured  by  2537  A  UV  absorption  at  the  shock  tube  test  station, 
assuming  an  absorption  coefficient  of  300/cm-atm.  From  the  difference 
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between  the  concentration  at  mixing  and  at  testing,  the  decomposition 
could  be  deduced.  These  measurements  gave  the  mole  fraction  of  03  ac¬ 
curate  to  approximately  0.05 $,  and  showed  decomposition  increasing  with 
time  from  0  to  30$  of  the  initial  03  concentration.  A  separate  study 
of  the  decomposition  process  by  IR  absorption  showed  that  at  dry  ice 
temperature  a  maximum  of  30$  of  the  03  lost  was  converted  to  C02,  and 
this  amount  of  C02  can  be  assumed  to  have  been  present  as  an  impurity 
in  the  test  gas  samples. 

Before  each  shock  tube  test,  the  tube  was  evacuated  and  then  filled 
with  a  passivation  sample  of  pure  03.  After  approximately  10  minutes, 
the  tube  was  again  evacuated  to  10“8  torr,  at  which  point  the  combined 
leak  and  outgassing  rate  was  about  2  x  10“4  torr/minute.  The  tube  was 
sealed,  filled  with  test  gas,  and  the  shot  fired  within  3  minutes.  As 
noted  above,  the  actual  03  concentration  in  the  test  gas  sample  was  mea¬ 
sured  by  UV  absorption,  and  total  sample  pressure  was  measured  both  by 
calibrated  Wallace  and  Tiernan  gages  and  by  a  sulfuric  acid  manometer. 
Total  pressures  varied  from  2  to  150  torr  and  were  measured  to  an  un¬ 
certainty  of  less  than  5$.  The  shock  velocity  was  measured  with  thin 
film  heat  transfer  gages  with  an  uncertainty  of  less  than  2$.  The  com¬ 
bined  uncertainties  of  the  physical  quantities  (pressure,  temperature, 
relaxation  time,  and  the  like)  were  approximately  ±  20$. 


Ill  OBSERVATIONS 


To  establish  baselines  for  the  effects  of  0  atoms  on  the  vibrational 
relaxation  of  N2,  a  series  of  preliminary  tests  were  made  on  Na-1$  CO, 

N2-5$  CO,  and  N2-l$C0-5#  Oa  mixtures.  Results  from  these  experiments 
are  presented  in  Figure  1,  together  with  curves  representing  earlier 
results  for  similar  mixtures  obtained  by  Millikan  and  White. ®>s  The 
"pure"  N2  results  can  be  seen  to  agree  very  well  over  their  common  tem¬ 
perature  range,  although  there  is  a  deviation  from  a  straight-line  ex¬ 
trapolation  at  temperatures  below  about  2000  K.  Such  a  deviation  could 
be  caused  by  small  concentrations  of  impurity  gases,  and,  in  part,  by 
the  CO  gas.  The  latter  effect  is  demonstrated  by  the  different  results 
for  mixtures  containing  1  and  5$  CO.  The  results  for  the  N2-02  mixture 
also  show  good  agreement  with  the  earlier  results.9  The  N2-02  measure¬ 
ments  of  this  study  were  limited  to  temperatures  below  2600  K  by  extra¬ 
neous  radiation  that  obliterated  the  relaxation  profile.  To  study  this 
problem,  some  experiments  also  were  conducted  with  radiation  filters 
that  passed  shorter  wavelengths;  in  these  tests,  the  extraneous  radiation 
occurred  at  even  lower  temperatures,  leading  to  speculation  that  the  radia¬ 
tion  is  from  C02  formed  by  reaction  between  CO  and  02.  Another  possibility 
is  that  the  radiation  arises  from  H20  formed  by  combustion  at  the  contact 
surface,  and  the  radiation  diffuses  forward  through  the  CO. 

In  the  N2-C0  and  N2-C0-02  tests,  the  CO  radiation  showed  a  single 
exponential  risetime  characteristic  of  an  IR  tracer  method  in  which  the 
vibrational  temperature  of  the  tracer  gas  is  completely  equilibrated  with 
that  of  the  N2.  When  03  is  added  to  the  test  gas  mixture,  however,  the 
radiation  exhibits  a  double  exponential  behavior  as  shown  in  a  typical 
oscilloscope  trace  and  a  semilogarithmic  plot  of  that  trace  in  Figure  2. 
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FIGURE  1 


VIBRATIONAL  RELAXATION  TiMES  FOR 
N-  -  CO  AND  N2  -  CO  -  0-  MIXTURES 
BY  THE  IR  TRACER  METHOD 


LABORATORY  TIME  ~  nsec 


FIGURE  2 
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TYPICAL  INTENSITY  HISTORY 
(MACH  4.8,  p,  =  45  torr,  T?  =  1530  K) 
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This  type  of  behavior,  as  observed  and  discussed  previously  by  Taylor, 
Camac,  and  Feinberg, 12  and  by  Sato,  Tsuchiya  and  Kuratani13  normally 
indicates  that  the  CO  has  relaxed  independently  to  a  considerable  ox- 
tent  before  it  couples  to  the  Na.  The  initial  risetime  Tj  then  gives 
the  rate  of  relaxation  of  CO  in  the  mixture.  As  discussed  in  a  follow¬ 
ing  section,  however,  this  does  not  seem  to  be  the  complete  explanation 
in  these  tests.  The  second  time  constant  Ta  in  nil  cases  is  related  to 
the  N2  relaxation  rate.  A  very  limited  number  of  tests  were  conducted 
with  an  Ar-1$  C0-1$03  mixture.  In  these  tests,  the  radiation  again 
showed  a  single  exponential  risetime  of  CO  relaxation  governed  by  the 
effects  of  0  atoms. 

A  series  of  tests  was  conducted  at  each  of  the  nominal  03  con¬ 
centrations  of  0.2,  0.5,  and  1.0$.  The  experimental  mixture  relaxation 
times,  pt2j  for  these  nominal  concentrations  are  presented  in  Figure  3. 
As  observed  also  by  Breshears  and  Bird,2  even  these  very  small  con¬ 
centrations  of  03  (and  consequently  0  atoms)  produce  substantial  reduc¬ 
tions  in  the  observed  Na  relaxation  times  (e.g.,  nearly  two  orders  of 
magnitude  for  1$  03  at  1000  K).  Because  it  was  not  always  possible  to 
add  exactly  the  desired  amount  of  03  to  the  test  gas  mixtures,  and 
because  of  decomposition,  the  actual  03  concentrations  deviated  from 
the  nominal  values  by  up  to  50$.  This  contributes  to  the  scatter  of 
the  data  in  Figure  3. 
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FIGURE  3  VIBRATIONAL  RELAXATION  TIMES  FOR 
n2  IN  A  MIXTURE  CONTAINING  OXYGEN 
ATOMS 
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IV  ANALYSIS 


Because  of  the  many  species  present,  the  number  of  vibration- 
translation  (V-T)  and  vibration-vibration  (V-V)  processes  taking  place 
in  the  gas  mixture  is  large.  We  can  summarize  the  possible  V-T  pro¬ 
cesses  as 


Na 

+  M  3S  Na*  +  M 

(7) 

°a 

* 

+  M  U  08  +  M 

(8) 

CO 

* 

+  M  5#  CO  +  M 

(9) 

C02 

+  M  si  C02*  +  U 

do) 

where  $  denotes  vibrational  energy  and  the  various  species  M  are  de¬ 
noted  in  the  following  discussion  as:  A  =  N2,  B  £  0a,  C  =  0,  D  =  CO, 
and  E  =  C02.  The  dominant  V-V  processes  are: 

♦  $  /X 

N,  +  00  58  Na  +  CO  (11) 

±  ± 

Na  +0aJiNa+0a  (12) 

♦  *  ,  . 

Na  +  C08  Na  +  CO,  .  (13) 

Because  of  their  small  relative  concentrations,  we  will  neglect  V-V 
interactions  between  CO,  02,  and  C02. 

The  equations  governing  the  vibrational  relaxation  of  gas  mixtures 
were  derived  by  Schwartz,  Slawsky,  and  Herzfeld, 14  and  they  subsequently 
were  put  in  a  more  convenient  form  by  Taylor,  Camac,  and  Feinberg. 12  In 
the  latter  form,  the  quantities  and  are  defined,  with  =  (E^-  E^)/ 
E^  where  E^  is  the  vibrational  energy  of  species  i  at  time  t,  is  the 


B-12 


50 


equilibrium  value,  and  is  the  mole  fraction  of  species  i.  The  per¬ 
tinent  equations  for  this  study,  namely,  for  the  relaxation  of  N2  and 
CO,  are: 


_S_  V1-^ 

v  l-exp(-©a 
T .  _  *—  i 


iAi-i  )  "I 

~ — . r  ■' 1  i  •  [l-exp(-©  )] 

l-exp(-0  )  L  A 

B  - 


v1-^)  w  ~i  r,  . 

l-exp(  -©^)  l-exp(-©D)  *  ^  exp  "a  ^ 


n 

V  1-1 
Taw  L 


«e(1-»a> 
-exp( -©. , 


Ui-O 


‘A* 

1-exp 


V  ~ 


[l-exp(-©A)]  (14) 


XA  oxp(©D)  [“  $D(1-$A)  $Ad-$D)  I 

+  v  exp(©  )  I  l-exp(-0  )  l-exp(-0  )  K1_exP(“®D)] . 

tad  A  A  D  J 

(15) 

In  these  equations,  ©  =  hv  /k,  t  is  the  V-T  relaxation  time  of  species 

^  V 

i  infinitely  dilute  in  species  j,  and  t  is  the  V-V  relaxation  time 

between  species  i  and  j.  The  first  set  of  terms  on  the  right  hand  side 

of  these  equations  thus  represents  V-T  relaxation,  while  subsequent  terms 

represent  V-V  transfer  with  each  of  the  possible  collision  partners. 

Note  that  in  the  case  of  V-V  transfer  between  C0a  and  N2,  we  restrict 

attention  to  the  v3  node  of  C02  (u)  =  2349  cm"1). 

6 
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Sato,  et  al.,  have  solved  a  system  of  equations  like  (14)  and  (15) 
when  there  was  only  V-V  exchange  between  two  species.13  In  principle,  a 
closed- form  solution  can  also  be  obtained  when  there  are  additional  V-V 
terms,  but  in  practice,  the  algebra  becomes  unduly  cumbersome.  We  will 
therefore  neglect  for  the  moment  the  V-V  contributions  of  both  02  and 
C02  and  will  show  later  that  these  terms  are  indeed  negligible.  We 

(V* 

also  assume  at  this  point  that  ©.  *  ©  and  lump  the  coefficients  of 

AD 

each  of  the  V-T  terms  into  a  single  parameter  T"1  ,  where 


X 

=  S  -i- 
j  Tij 


(16) 


The  equations  then  simplify  to 


(14a) 


and 


d$ 


dt 


—  $.  + 
v  A 

TAD 


(15a) 


The  general  solution  of  these  equations  was  given  by  Sato,  et  al.,  *8 

their  equations  (9)  through  (l3),13  and  it  need  not  be  repeated  here. 

When  X  »  X  and  T  <  T  ,  the  solutions  simplify  considerably  and  as 
A  B  DA 

given  by  Sato,  et  al.,  are: 


§A  =  exp(-t/TA) 


(17) 
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**JUr-*‘ 


and 


§  = -  exp[-( -  -  —-)t]  +  — —  exp(-t/T  )  , 

D  1+r  XT  T.  J  1+r  '  A'  * 

D  A 


(18) 


where 


r  =  l/C(T-  -  V‘)  W 


(19) 


The  solution  f-  5.  indicates  that  N„  will  relax  with  a  single  exponential 
A 

risetime  governed  only  by  V-T  processes,  but  that  for  shows  that  CO 
will  show  a  double  exponential  behavior  as  discussed  in  the  previous  sec¬ 
tion.  When  «  Td,  r  is  large  and  the  first  term  of  equation  (18) 
becomes  negligible.  This  is  normally  the  case  in  the  IR  tracer  method. 

Examination  of  equation  (18)  indicates  that  the  initial  time  con¬ 
stant  Tj  observed  in  the  relaxation  profile  is  T^,  and  that  the  final 
time  constant  Ta  is  T  .  Furthermore,  the  parameter  r,  which  can  be 
determined  by  extrapolating  the  final  relaxation  profile  back  to  time 
zero,  can  be  used  in  equation  ( 19)  to  determine  T  once  T  is  known, 

U  A 

since  t._  has  already  been  determined  by  several  investigators. i3> 15> 18 

Ai/ 

Since  T  can  thus  be  determined  in  two  ways,  we  have  a  means  of  check- 
D 

ing  on  the  internal  consistency  of  the  present  results  for 

When  Td  «  T^,  the  simplified  solution  above  may  not  adequately 
account  for  V-V  transfer  from  CO  to  N2.  In  this  case,  the  solution  of 
Sato,  et  al.,  can  be  modified  to  give,  instead  of  equations  (17)  and 
( 18),  the  expressions: 
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Since  (r/l+r)  <  1  and  (X  /x  )  =  0.01,  we  see  from  equation  (20)  that 

D  A 

V-V  transfer  from  CO  to  N2  becomes  important  when  T  ^0.1  T  .  We  find 

JJ  A 

from  measurements  of  pt  _  and  from  calculations  that  the  differences 

CO-O 

between  T  and  T  are  in  some  cases  greater  than  this,  and,  therefore, 
DA 

we  use  these  more  complete  equations  to  calculate  the  rate  of  relaxation 
of  N2  in  this  study. 

Equation  (20)  can  also  be  used  to  examine  the  effects  of  02  and  C02 
on  the  relaxation  of  Na.  We  have  already  made  the  assumption  that  each 
of  these  species  relaxes  by  V-T  processes  or  by  V-V  interaction  only  with 
N2.  In  each  case,  then,  the  effect  of  species  i  on  N2  will  be  negligible 
if 


x<  Ta  Ta 

n  -  1  A  r  ,  i  A 

i  "  XA  T  1+r  XA  m  ,  v 
A  i  A  T1+taj 


«  1 


(22) 


We  can  calculate  T  ,  including  the  effects  of  0  atoms,  from  the  results 

,  v 

of  Kiefer  and  Lutz, 1  and  t  from  the  results  of  White9  and  Kamimoto 

Na-  o2 

and  Matsui. 10  In  this  case,  the  V-V  transfer  is  relatively  slow  (r  «  l), 
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and  since  X  £  0.015,  we  find  R  <  0.05.  For  CO,  ,  preliminary  informa- 
°2 

tion  on  the  relaxation  of  CO-  by  0  atoms  indicates  that  process  to  be 


negligible  here.17  The  relaxation  times  T  and  t 


CO- 


Na-COa 


can  then  be  ob¬ 


tained,  for  example,  from  the  results  of  Taylor  and  Bitterman. 18  Taking 

the  maximum  concentration  of  CO-  as  X _  £  0.001  then  gives  R_  <  0.08. 

a  C0a  C0a 

Therefore,  within  the  accuracy  inherent  in  this  study,  the  effects  of  0a 
and  C02  can  safely  be  neglected. 


Having  related  Tx  and  t2  to  T  and  T  through  equation  (21),  we  can 

U  A 


proceed  to  find  t„  „  and  t _  ^ 

Na-  0  CO-O 


From  equation  (16),  we  have 


1 _ Nj, 


Na-Na 


X. 


'n2-o2 


Na-0 


CO 


CO, 


Na-C0 


Na-COa 


(23) 


and 


N„ 


CO-N, 


X. 


'JL. 


CO-N, 


X. 


CO-O 


CO 


CO-CO 


CO, 


co-coa 


(24) 


However,  on  the  basis  of  Millikan  and  White's  correlation, 6  for  example, 

we  assume  t  t  and  t  __  approximately  equal  to  t„  „  , 

and  likewise  t  „  ,  T  _  ,  and  t  __  equal  to  t  „.  Then,  since 
CO-Na  CO-Oa  C0-C02  CO-CO 

Xq  £  0.01  so  that  XN  +  X^  +  Xc()  +  X(X)  =  1,  we  obtain  the  expressions 


Na-  Na 


Na-  0 


N.-  Na 


-  T. 


(25) 
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and 


TC0-0 


X  •  T  *  T 

0  CO-CO  D 

T  -  T 
CO-CO  D 


(26) 


The  times  t  .  „  are  taken  from  the  curve  of  Figure  1, 

.2™  ^*2 

are  taken  from  the  results  of  Hooker  and  Millikan. 19 


and  the  x _  's 

CO-CO 


In  analyzing  the  data,  thermodynamic  properties  behind  the  shock- 
wave  were  calculated  for  the  gas  mixture,  pressure,  and  Mach  number'of 
each  test.  The  fraction  of  03  converted  to  0  atoms  was  also  calculated, 
using  the  method  described  in  the  Appendix.  Results  for  each  test  are 
presented  for  a  temperature  that  is  an  average  between  the  values  cal¬ 
culated  for  vibration  frozen  and  for  full  vibrational  equilibrium,  with 
the  03  assumed  dissociated  in  both  cases.  The  maximum  difference  be¬ 
tween  these  two  temperatures  is  approximately  12$. 
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V  RESULTS 


In  view  of  the  possible  contribution  of  V-V  transfer  from  CO  to  the 

relaxation  of  Na  as  indicated  in  equations  (20)  and  (21),  we  report  first 

our  measured  relaxation  times  pr  _  •  The  results  derived  from  the 

CO-O 

Ar-1$  C0-1$  03  measurements,  presented  in  Figure  4,  indicate  that  PT£q_q 
is  approximately  0.04  atm-|is.  Although  limited  in  number,  these  values 
are  in  good  agreement  with  the  more  extensive  measurements  of  Center,  which 
were  obtained  by  using  similar  gas  mixtures.17  Unfortunately,  the  re¬ 
sults  for  pt^  Q  obtained  from  the  Na-C0-0s  tests  are  neither  consistent 
nor  readily  explainable.  As  discussed  above,  pt^  Q  can  be  determined 
both  from  the  initial  risetime  and  from  the  intercept  parameter  r. 

The  results  based  on  the  initial  risetime  tq  are  on  the  order  of  ten 
times  longer  than  those  presented  in  Figure  4,  while  those  based  on  the 
parameter  r  range  from  l/lO  to  l/4  as  long  as  those  of  Figure  4.  From 
these  anomalies,  the  conclusion  is  that  neither  the  initial  risetime  nor 
the  intercept  are  solely  related  to  the  CO-O  relaxation  process;  nonce 
for  the  time  being,  this  interesting  effect  must  await  an  alternative 
explanation. 


Because  of  the  uncertain  interpretation  of  the  initial  risetime, 
we  calculate  the  relaxation  time  of  CO  in  the  Na-C0-03  mixture  using  the 
results  for  from  Figure  4.  The  relaxation  times  p>rN  for  Na  by 

0  atoms  are  then  determined  from  equations  (21)  and  (25),  and  the  results 


are  presented  in  Figure  5  for  each  of  the  nominal  03  concentrations.  In¬ 
cluding  the  uncertainties  of  the  measurements  of  the  physical  quantities 
(±  20$),  of  the  correction  for  V-V  transfer  from  CO  and  of  possible  con¬ 


tributions  from  other  minor  species  or  impurities,  the  total  uncertainty 
of  these  results  is  estimated  to  be  ±  50$ 


Also  shown  in  Figure  5  are  the  earlier  results  of  Breshears  and 
Bird,2  the  low  temperature  results  of  McNeal,  et  al.,5  and  the  predictions 
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FIGURE  5  RELAXATION  TIMES  FOR  THE  VIBRATIONAL  RELAXATION  OF  N,  BY  OXYGEN  ATOMS 


of  Fisher  and  Bauer.7  Although  the  data  show  substantial  scatter,  the 
current  results  agree  very  well  in  magnitude  with  both  the  earlier  shock 
tube  results  and  with  the  low-temperature  results. 

As  noted  earlier,  a  very  large  disagreement  exists  between  the 
temperature-dependences  of  the  theoretical  and  experimental  results.  In 
the  theoretical  prediction,  the  temperature-dependence  results  from  a 
specification  of  the  height  of  the  crossing  point  of  the  N2  +  0(3P)  and 
the  N2  +  0(1D)  potential  curves  above  the  asymptotic  limit  of  the 
N2  4-  0(3P)  curve.  Fisher  and  Bauer  used  a  value  of  0.8  eV  for  this 
barrier  height  or  "activation  energy"  based  on  the  unimolecular  decom¬ 
position  of  N20. 7  If  their  model  is  to  apply,  the  activation  energy 
must  actually  be  less  than  0.2  e\J  although  the  theoretical  curve  does 
not  fit  the  shape  of  the  experimental  data  very  well  even  under  that 
assumption.  The  activation  energy  could  differ  from  the  deduced  value 
of  0.8  eV  because  of  vibrational  excitation  effects  in  the  N20  dissoci¬ 
ation  or  because  of  orientation  effects  in  the  0-Na  interaction.7  It  is 
clear  that  additional  work  is  required  to  provide  an  adequate  theoretical 
explanation  of  these  rapid  relaxation  rates. 
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VI  CONCLUSIONS 


The  results  presented  here  for  the  rate  of  relaxation  of  N2  by  0 
atoms  are  important  in  that  they  serve  to  substantiate  both  the 
existing  high-  and  low-temperature  results.  In  particular,  with  this 
verification  of  the  low-temperature  results,  the  observed  vibrational 
temperature  of  N2  in  the  upper  atmosphere  can  now  be  explained. 20 > 21 
The  pronounced  difference  between  the  observed  and  predicted  temperature 
dependence  of  this  relaxation  rate  also  presents  a  challenge  to  find 
an  adequate  theoretical  understanding  of  this  process. 

The  relaxation  rate  for  CO  by  0  atoms  determined  from  tests  in 
Ar-C0-03  mixtures  is  in  good  agreement  with  other  recent  measurements. 
This  rate  is  approximately  25  times  as  fast  as  the  N2-0  rate  and,  co¬ 
incidentally,  is  approximately  the  same  as  that  for  the  relaxation  of 
02  by  0  atoms.  The  greater  rapidity  of  the  CO-O  and  02-0  relaxation, 
compared  with  that  of  Na-0,  may  be  the  result  of  atom-exchange  processes 
in  the  former  cases. 
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Appendix 


DECOMPOSITION  OF  SHOCK-HEATED  OZONE 

The  decomposition  of  shock-heated  0..  and  its  conversion  to  0  and 
02  is  governed  by  the  reactions: 

03  +  M  ^  0  +  02  +  M  (A-l) 

and 

0  +  03  -2  2  02  .  (A-2) 


Wray22  has  shown  that,  at  the  completion  of  the  decomposition  process, 

the  mole  fraction  of  0  atoms  (X  )  and  of  initial  0,(x  )  are  related  b*' 

0  3  03 

X0  =  (1H5)  xo  '  (A"3) 

3 

where  g  is  given  by  the  solution  of  the  equation 

in  [l-G(l-g)]  +  |  G(l-g)  +  \  G  =  0  .  (A-4) 


The  parameter  G  is  related  to  X  and  the  reaction  rates  k2  and  k2  by 


G  =  X  (kjkj  . 


(A-5) 


It  can  be  shown  that  g  is  asymptotic  to  G  for  G  «  1  and  to  1-l/G  for 
G  »  1.  With  this  information  and  some  numerical  solutions  of  equation 
(A-4),  the  variation  of  g  with  G  for  all  values  oi'  G  has  been  deter¬ 
mined  and  is  presented  in  Figure  A-l.  As  noted  earlier,  for  high  tem¬ 
peratures  where  k2/kx  =  1  and  for  X  «  1,  g  -»  G  -*  0  and  there  is  a 
one-to-one  conversion  of  03  to  0  atoms.  Conversely,  for  low  temperatures 
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SA-1 664-7 

FIGURE  A-1  THE  OZONE  DISSOCIATION  PARAMETER 

0  VERSUS  THE  PARAMETERS  G  AND  1  -  1/G 

(8  -  1  -  Xn/Xn  ,  G  *  Xn  •  k,/k,) 


where  k2/kt  »  1,  G  may  be  much  larger  than  1  so  that  (3  -*  1.  Then  there 
are  very  few  0  atoms  produced. 

The  decomposition  of  03  has  been  studied  by  Jones  and  Davidson33 

2  6 

(JD),  Myers  and  Barlle24  (MB),  Benson  and  Axworthy25  (BA),  and  Michael. 
They  have  given  the  following  rate  equations,  in  units  of  liter/mole-s: 

( JDl)  :  kj  =  5.8X1011  exp(-23,500/RT) 

1  1,75 

( JDIl)  :  kt  =  3.4X108T^(24,350/RT)  exp(-24350/RT) 

(JD)  :  k2  =  2.4X1010  exp( -5,600/RT) 

(MB)  :  k2  =  1. 04Xl09T^exp( -5, 288/RT) 

(BA)  :  k2  =  3.37xlOi0  exp(-5,700/RT). 

The  magnitudes  of  these  rates  are  such  that  the  largest  ratio  of  k2/kt 
is  formed  from  the  combinationka(MB)/k1(jDIl)  while  the  smallest  ratio 
is  formed  from  k2(JD)/k1(jDl),  as  shown  in  Figure  A-2.  The  latter 
ratio,  based  on  data  in  the  680-920  K  temperature  range,  has  been  the 
standard  of  comparison  for  most  studies  involving  03  decomposition-in 
the  shock  tube,  and  most  investigators  have  found  that  k2/k1  exceeds  the 
JD  extrapolation  by  an  increasing  amount  as  the  temperature  increases. 

For  example,  Michael  reports  a  factor  of  2.5  at  1300  K,26  Wray  estimates 
a  factor  of  10  for  1300  to  2900  K, 32  and  Kiefer  and  Lutz  estimate 
values  of  3  to  30  for  1500  to  3000  K. J^27  The  combination  k3(Jffl)/k1(JDIl) 
comes  close  to  giving  the  observed  ratio  for  all  temperatures,  but  there 
is  little  theoretical  justification  for  those  rates.  Rather,  it  is  more 
likely  that  k2  is  well  represented  by  the  JD  or  BA  expressions  while  kj 
deviates  from  all  of  the  given  expressions  as  temperature  increases.  In 
this  work,  therefore,  we  have  chosen  to  use  the  k2(BA)  rate,  together 
with  kt(jDl)  modified  downward  by  a  factor  (900/t)2,  which  gives  agree¬ 
ment  with  the  empirical  observations.  The  ratio  resulting  from  this 
combination  is  also  given  in  Figure  A-2. 
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The  fractional  decomposition  of  0g  into  0  atoms  determined  from 
equation  (A-4)  and  the  various  ratios  kj/kj  of  Figure  A-2  is  given  in 
Figure  A-3  for  the  cases  X  =  0.002  and  0.01.  In  the  case  of  X  =  .002, 

°3  °3 

the  differences  resulting  from  the  extreme  choices  of  ka/k1  range  from 
20$  at  1000  K  to  5$  for  T  >  2000,  while  for  X =  .01  the  differences 

°3 

range  from  45$  at  1000  K  to  15$  at  2000  K.  The  results  using  the 
empirical  ratio  of  kjj/kj  are  regarded  as  a  best  estimate  of  Xq/Xq  ,  and 
the  uncertainty  is  probably  somewhat  less  than  the  maximum  difference 
obtained  using  the  limiting  ratios  of  k2/kj. 
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